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Abstract 

The NSl protein of influenza A virus (lAV) is a multifunctional virulence factor. We have previously characterized gain-of- 
function mutations in the NSl protein arising from the experimental adaptation of the human isolate A/Hong Kong/1/ 
1968(H3N2) (HK) to the mouse. The majority of these mouse adapted NSl mutations were demonstrated to increase 
virulence, viral fitness, and interferon antagonism, but differ in binding to the post-transcriptional processing factor 
cleavage and polyadenylation specificity factor 30 (CPSF30). Because nuclear trafficking is a major genetic determinant of 
influenza virus host adaptation, we assessed subcellular localization and host gene expression of NSl adaptive mutations. 
Recombinant HK viruses with adaptive mutations in the NSl gene were assessed for NSl protein subcellular localization in 
mouse and human cells using confocal microscopy and cellular fractionation. In human cells the HK wild-type (HK-wt) virus 
NSl protein partitioned equivalently between the cytoplasm and nucleus but was defective in cytoplasmic localization in 
mouse cells. Several adaptive mutations increased the proportion of NSl in the cytoplasm of mouse cells with the greatest 
effects for mutations Ml 061 and D125G. The host gene expression profile of the adaptive mutants was determined by 
microarray analysis of infected mouse cells to show either high or low extents of host-gene regulation (HGR or LGR) 
phenotypes. While host genes were predominantly down regulated for the HGR group of mutants (D2N, V23A, F103L, 
M106I-FL98S, L98S, M106V, and M106V-FM124I), the LGR phenotype mutants (D125G, IV1106I, V180A, V226I, and R227K) were 
characterized by a predominant up regulation of host genes. CPSF30 binding affinity of NSl mutants did not predict effects 
on host gene expression. To our knowledge this is the first report of roles of adaptive NSl mutations that impact 
intracellular localization and regulation of host gene expression. 
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Introduction 

The influenza A virus (lAV) NSl protein possesses multiple 
fijnctions that support virus replication. NS 1 can engage in many 
functions due to its ability to translocate to both the nucleus and 
cytoplasm of infected cells and interact with numerous cellular and 
viral factors including RNA. Cytoplasmic activities include 
blocking viral RNA detection by RIG-I signalling of type I 
interferon (IFN) induction as well as inhibition of IFN-stimulated 
antiviral proteins, suppression of the host cell apoptotic response, 
and enhancement of viral protein synthesis [1]. In the nucleus, 
NS 1 binds cellular post-transcriptional processing factors including 
the cellular post-transcriptional host factor cleavage and polyad- 
enylation specificity factor 30 (CPSF30), which has been reported 
to result in a blockade of host gene expression, including type I 
IFN [2,3]. 



NSl localization is governed by two known nuclear localization 
signals (NLS) (NLSl: aa 34-38, NLS2: aa 203-237), which 
interact with the cellular protein importin a. [4] and induce rapid 
nuclear localization following translation [5]. Later in infection, 
the protein is detected in both the nucleus and the cytoplasm [6], 
which is attributed to interaction of its NES element (137-147) [6] 
with the nuclear pore complex. Cytoplasmic NSl levels are also 
influenced by an inhibitory sequence (148-161) adjacent to the 
NES [7] . NS 1 also localizes to the nucleolus and contains a NoLS 
involving key basic residues Arg-224 and Arg-229 [4,8], however 
the role of NSl nucleolar localization and its contribution to viral 
replication is unknown. 

The ability of the NSl protein to translocate within discrete 
cellular compartments and engage in pro-viral functions including 
blockade of host gene expression and type I interferon (IFN) 
induction has been shown to be essential for optimum virus 
replication. Alanine substitutions at NLS residues 38 and 41 were 
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observed to reduce NS 1 nuclear localization and the ability to bind 
RNA, resulting in increased IFNot/p production and attenuation 
of virus replication [9,10]. Li et al. (2011) showed alanine 
substitutions at Leu-69 and Leu-77 within the highly conserved 
NSl linker region induced de novo localization to the nucleolus, and 
dramatically reduced the ability of NSI to limit IFN induction 
in vitro [11]. In addition, NSl binding affinity for nuclear and 
nucleolar factors, including CPSF30, nucleolin, and fibrillarin, also 
influence protein's sequestration in the nucleus and nucleolus, 
respectively [8,12]. In the case of CPSF30, amino acid mutations 
F103S and M1()6I disrupted CPSF30 binding while shifting NSl 
localization to the cytoplasm [12]. 

Although intracellular localization of the NSl protein is well 
documented, the effect of virus adaptation to a novel host on NSl 
protein localization remains unknown. Cellular distribution of the 
NSl protein is dependent upon virus-host interactions; several 
reports have shown strain-specific and cell-line specific localization 
patterns of the NSl protein [4,12,13]. Moreover, the nuclear 
import signals of influenza RNP (polymerase subunits (PBl, PB2, 
and PA) and NP) have been shown to be determinants of adaptive 
evolution in mice [14—17] and for highly pathogenic avian 
influenza in mammalian cells (reviewed in [18]). We hypothesize 
that subcellular trafficking of NSl is a determinant of host 
adaptation. 

The influenza A virus host range is extensive; while aquatic 
waterfowl species define the primary viral reservoir, many avian 
and mammalian species are also hosts. However, the genetic 
determinants of lAV host range and virulence remain largely 
unknown. To this end, we previously reported the experimental 
evolution of the human influenza A virus isolate A/Hong Kong/ 
1/1968 (H3N2) (HK) to high virulence in the mouse [15,17] and 
identified multiple mouse-adapted (MA) mutations in the PB2 
[16,17], HA [16,17,19] as well as the NSl protein [20,21]. These 
mutations, including 11 in the NSl gene, were confirmed to 
mediate host adaptation, assessed by enhanced virus replication in 
mouse lungs and/or mouse cells [21]. These mutants possessed 
additional adaptive phenotypes, including enhanced viral gene 
expression (at both transcriptional and translational levels) as well 
as an ability to reduce IFN-P induction in mouse lungs. 
Rc-markably, 8 mutants lost the ability to bind CPSF30 while 4 
maintained binding [21], and this trait did not correlate with IFN- 
P production in vivo in the mouse lung. One of the NSl mutants, 
D125G, also induced expression of a novel viral protein, the non- 
structural 3 protein (NS3), due to activation of an alternative splice 
site [21,22]. Although the function of this novel protein is currentiy 
under investigation, we have confirmed that expression of the NS3 
protein enhances viral replication for the recombinant HK virus 
expressing the NSl D125G mutation [22]. 

The adaptive phenotypes associated with MA NSl mutants 
were thus wide ranging, but more importantly, involved NSl 
protein functions occurring in both the nucleus and the cytoplasm 
of infected cells. As the HK virus is a human virus isolate, we 
hypothesized that its evolution to high virulence in the mouse 
selected for mutations in the NSl protein that would favour NSl 
cellular localization optimal for replicative fitness in the mouse. 
Furthermore several human HK NSl mutations selected on 
mouse adaptation mapped to regions within or near to known 
NLS and NES domains. We next sought to determine the effect of 
MA NSl mutations on NSl localization in mouse and human 
infected cells, and whether altered NSl localization correlated 
with lAV adaptation to a novel host. To further understand the 
role of CPSF30 binding during infection, we also assessed whether 
CPSF30 binding affinity of the adaptive mutants correlated with 
altered host gene expression in mouse cells. 



We found that the MA NSl mutant proteins were present at a 
signrficantiy higher proportion in the cytoplasm of infected mouse 
cells compared to wild-type HK NSl, and that NSl mutants 
altered NS 1 cellular localization in a host cell-dependent manner. 
We previously have shown that the majority of MA NS 1 mutations 
abrogated binding to the CPSF30 protein [21]; here we 
demonstrate two distinct phenotypes of altered host gene 
expression with either high or low effects upon infection of mouse 
cells with NSl mutant HK viruses, and moreover, that CPSF30 
binding affinity did not correlate with reduced expression of host 
genes. To our knowledge, this study is the first to identify the roles 
for NSl localization in host adaptation as well as in regulating 
general host gene expression. 

Materials and Methods 

Cell Lines 

MDCK cells (Madin-Darby canine kidney; Health Canada), 
Ml cells (Mouse kidney epithelium; ATCC), and A549 cells 
(human lung carcinoma; ATCC) were maintained in minimum 
essential medium (MEM). MEM was supplemented with L- 
glutamine (2 mM), Penicillin (100 U/ml), Streptomycin (100 ug/ 
ml) and fetal bovine serum (FBS) (10%). All cells were incubated at 
37°C in the presence of 5% CO2. 

Viruses 

The H3N2 human influenza isolate A/Hong Kong/ 1/1968 
(HK-wt) was originally obtained from the Laboratory Center for 
Disease Control (Health Canada, Ottawa, Canada). The HlNl 
human influenza isolate A/Puerto Rico/8/34 (PR8) was obtained 
by reverse genetics of the 8 viral genome segments expressed on 
pHW2000 plasmids gratefully obtained from R. Webby (St. Jude's 
Children's Research Hospital, Memphis). Th(; laboratory adapted 
HlNl influenza isolate A/WSN/33 (WSN) was obtained by 
reverse genetics of the 8 viral genome segments in pHH21 
plasmids obtained from Yoshihiro Kawaoka that were subcloned 
into the pLLB plasmids for virus rescue (University of Madison, 
Wisconsin). Viruses were plaque purified on MDCK cells and 
grown in the allantoic cavity of 10-day-old specific pathogen free 
(SPF) chicken embryos (Canadian Food Inspection Agency, 
Ottawa). 

Mouse Adaptation of A/Hong Kong/1/1968 (H3N2) 

The generation of mouse-adapted (MA) variants of A/Hong 
Kong/ 1/1968 (H3N2) was described previously [15,17]. Briefly, 
highly virulent variants were selected after 20 or 21 passages 
(MA20 or MA21) in the CD-I mouse lung. Following fmal 
passage, lung tissue was homogenized then serially plaque purified 
on MDCK cells. 

Generation of Recombinant NS Mutant Hong Kong 
viruses 

To determine the role of each mutant MA NSl gene in 
virulence, recombinant HK viruses were generated to express the 
mutations of each mouse-adapted NS gene segment (rHK NS 
mutants) as previously described [21]. The NS gene segments of all 
rescued viruses were sequenced to verify that they were free of 
unwanted mutations. 

Confocal Microscopy 

Localization pattern of NS 1 in infected cells were observed by 
seeding Ml or A549 cells to 60-80% confluency on 8 chamber 
immunoslides then infecting the cells with rHK NSl mutants or 
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HK-wt virus, as well as the prototype mouse-adapted viruses PR8 
and WSN, at an MOI = 3 in the absence of trypsin. At 16 hpi, the 
cells were washed with 1 x PBS (PBS), permeabilized with 0. 1 % 
Triton XI 00 in PBS for 5 minutes, then washed with PBS and 
fixed in 3.7% formaldehyde for 10 minutes. Following two washes 
in PBS, the cells were incubated with rabbit anti-HK NSl 
polyclonal antibody, rabbit anti-HK NP polyclonal antibody, or 
mouse anti-NSl monoclonal antibody overnight at 4°C. AU 
antibodies were previously pre-adsorbed in mouse lung extract to 
reduce non-specific cross-reactivity and suspended in PBS washing 
buffer containing 3% bovine serum albumin and 0.3% Triton X- 
100. The cells were washed twice with PBS washing buffer then 
incubated with secondary Cy3-conjugated donkey anti-rabbit or 
donkey anti-mouse IgG, respectively (Jackson ImmunoResearch 
laboratories Inc., ME), for 2 hours at room temperature. Following 
two washes with PBS, the nuclei were labeled by DAPI incubation 
in the dark for 15 minutes (300 ng/mL). Excess stain was removed 
by 3 washes with dH20. Cover slides were mounted with Dako 
mounting medium (Dako, Denmark) as per the manufacturer's 
protocol. The cells were viewed using a Zeiss LSM 510 MET A/ 
AxioVert 200 Confocal Microscope 63 x oil immersion objective 
lens. Representative images were processed in parallel (GNU 
Image Manipulation Program 2.0). 

Quantification of NSl in Cellular Compartments of 
Infected Cells 

Confluent monolayers of Ml or A549 cells were infected (n = 2- 
3) with rHK NSl mutant or HK-wt viruses as well as the 
prototype MA viruses A/PR/8/34 and A/WSN/33 (MOI = 2) in 
the absence of trypsin. At 16 hpi, the cells were washed once with 
ice cold 1 X PBS then lysed with ice cold 1% NP-40 (in PBS) and 
scraped to collect adherent cells. Cells were fractioned as 
previously described [23,24]. Briefly, an aliquot of the cell lysate 
was obtained to represent the whole ceU lysate, and the remainder 
was centrifuged for 2 minutes at 13.5 K RPM. The supernatant 
containing the ceU cytoplasm was then coUected. The nuclear 
pellet (that employed twice as many ceUs as the W and C samples 
to counteract losses in preparation) was further washed with NP- 
40 lysis buffer and re-centrifuged. The supernatant was discarded 
and the pellet was re-suspended in 50 uL NP-40 lysis buffer. 
Whole cell, cytoplasmic, and nuclear fractions were then mixed 
with 4x Laemmli sample buffer and denatured by 3 minutes at 
90°C prior to protein separation by SDS-PAGE gel electropho- 
resis. Immunoblots were probed for NSl (polyclonal rabbit anti- 
NSl), histone H3 (rabbit anti-histone H3, CT, pan clone A35; 
Millipore, Billerica, MA, USA) and P-tubulin (mouse anti-tubulin 
a; Sigma, Burlington, ON). Antibodies were detected by 
incubation with HRP conjugated goat anti-rabbit or goat anti- 
mouse antibody (Sigma, Burlington, ON) foUowed by SuperSignal 
West Pico chemiluminescent substrate (Pierce, Thermo Fisher 
Scientific, Nepean, ON) and exposure to film. 

Densitometry Analysis 

Band density of proteins detected by Western blot was 
quantified using the UN-SCAN-IT software (Scientific Software 
Solutions). As the NSl bands associated with HK-wt and 
M106V+M124I infections were much \\x'akc'r than the remainder 
of NSl bands, relative NSl band densit)' was determined using 
two exposures (a longer exposure to quantify HK-wt and Ml 06V 
+M124I NSl bands, shown in each figure, and a shorter exposure 
that allowed the remainder of NSl bands to remain below 
oversaturation (data not shown)). To determine the cellular 
distribution of the NS 1 protein in infected ceUs, we determined 
the ratio of cytoplasmic NSl to nuclear NSl. As the histone H3 



band was inconsistent between the whole ceU and nuclear fractions 
(possibly due to losses from the washing of nuclei by sedimentation 
involved in the fractionation procedure), the amount of NSl 
protein in the nuclear and cytoplasmic fractions of a cell infected 
with a given virus was normalized to histone H3 and tubulin levels 
of the whole ceU lysate, respectively. The fraction of cytoplasmic to 
nuclear NSl protein was calculated by formula; (histone W/ 
histoneN) xNS 1 N = normalized NS 1 N to normalize nuclear NS 1 
quantities to W fraction, as weU as for cytoplasmic fractions 
(tubulin W/tubulin C)xNSl C = normalized NSl C, to obtain 
values to determine the ratio of C/N (using both values as 
normalized to W). 

Microarray Analysis 

Confluent monolayers of Ml cells in 35 mm dishes were 
infected with HK-wt or rHK NS mutant viruses (MOI = 2) or 
mock-infected with PBS in triplicate for n = 3 biological replicates 
as per the previous infection protocols. At 8 hpi, the cellular RNA 
was coUected using manufacturer's protocol (QIAGEN RNeasy 
mini-kit) quantified using a PowerWave XS2, Microplate Spec- 
trophotometer (BioTek Instruments, Inc. Winooski, VT United 
States ), and ribosomal RNA quality and integrity was confirmed 
using an Agilent 2100 Bioanalyzer (Agilent Technologies Canada 
Inc. Mississauga, Canada). The microarray analysis was per- 
formed by StemCore Laboratories (Ottawa, Ontario, Canada) 
using Genechip Mouse Exon 1.0 ST Array (Aflymetrix, Santa 
Clara, CA, USA). Microarray gene expression data were analyzed 
by Flexarrayl.6. 1 after normalization using Affymetrix Power 
Tools (AFT) [25]. Genes were considered as up or down regulated 
relative to mock infected cells if they were £1 or & 1 log2 fold 
difierent (^ or S2 fold differences) expression level and were 
statistically significant at the pS0.05 level using ANOVA. Genes 
were included for hierarchical clustering analysis among mutants if 
they were differentially regulated (S or >2 fold differences) and 
significandy different from mock-infected cells for one or more of 
the mutants. Microarray data were further analyzed using 
Bioconductor software packages in R for clustered heat maps of 
gene expression and mutants, Pheatmaps, (pheatmap_0.7.4.tar.gz) 
analysis using the complete linkage and Euclidian distance 
method; and gene ontology', GoProfiles (goProfiles_1.20.0.tar.gz), 
(avaflable at http://bioconductor.0rg/biocLite.R) with outputs 
plotted with gplots_2.1 l.O.tgz and user generated software. The 
microarray data have been deposited in NCBI's Gene Expression 
Omnibus and are accessible through GEO Series accession 
number GSE48217 (http://www.ncbi.nlm.nih.gov/geo/query/ 
acc.cgi?acc = GSE48217). 

Statistical Analysis 

Unless otherwise indicated, statistical significance was measured 

using the student t-test (Microsoft Excel, 2007) using the 
parameters of equal variance and two-tailed significance with a 
probabilit}' of <0.05 considered as statistically significant. Values 
were calculated as means ± standard deviation (SD), unless 
otherwise noted as means ± standard error (SE). 

Results 

NSl Mutants Alter Subcellular Localization of the NSl 
Protein in Infected Mouse Cells 

We have previously shown that NSl mutations selected upon 
mouse adaptation of the A/Hong Kong/1/1968 H3N2 virus (MA 
NSl mutations) are adaptive, and associated with increased 
virulence in the mouse, decreased IFN|3 production in the mouse 
lung, and increased viral repUcation, both in vivo and in vitro, of 



PLOS ONE I www.plosone.org 



3 



December 2013 | Volume 8 | Issue 12 | e84673 



Adaptation and Cytoplasmic Localization of NS1 



which the latter was associated with enhanced viral gene 
expression at the level of viral mRNA and protein synthesis 
[21]. We next assessed whether virus adaptation to a novel host 
alters NS 1 localization using confocal microscopy. Mouse M 1 cells 
were infected with recombinant HK viruses expressing MA NSl 
mutations (rHK NSl mutants) at a MOI of ,3 foUowxxi by 
formaldehyde fixation at 16 hours post infection (hpi) and 
immunostaining using an anti-NSl or an anti-NP (for virus 
protein control) polyclonal antibody with a Cy3- conjugated 
secondary antibody, and DAPI staining to detect the nucleus with 
images obtained by confocal microscopy using a 63 x oil 
immersion objective. Given that mouse adapted NSl mutants 
were previously shown to differentially increase viral fitness and 
viral protein production in the mouse [21], we saw similar or 
increased intensity of NP staining for each mutant relative to HK- 
wt (Figure SI). 

The majority of mouse cells that were infected with the HK-wt 
virus demonstrated staining for the NS 1 protein exclusively in the 
nucleus (Figure 1), with only 6% of cells positive for NSl in both 
cellular compartments (Figure S2). However, mouse cells infected 
with NSl mutants D125G and M106I+L98S, and to a lesser 
extent M106I and V226I, significantly increased the observed 
frequency of NSl staining within the cytoplasm to up to 64% 
(p<0.001; D125G) (Figure 1 and S2). Thus several NSl adaptive 
mutations increased cytoplasmic localization of the NS 1 protein in 
mouse cells detectable by microscopy. 

NSl Mutants Temporally Increase Localization of the NSl 
Protein to the Cytoplasm in Infected Mouse Cells 

Because several MA NSl mutants altered NSl protein 
locahzation in mouse cells at 16 hpi, we next assessed the pattern 
of cytoplasmic staining at 4, 8, 16 and 24 hpi to determine the 
temporal pattern of NSl subcellular localization. We analyzed the 
NSl mutants D125G and M106I+L98S, which exhibited the most 
dramatic cytoplasmic staining phenotype, as well as two prototype 
laboratory LAV strains that had been previously adapted to the 
mouse, A/WSN/1933 (HlNl) (WSN) and A/PR/8/1934 (HlNl) 
(PR8). Representative confocal images for each time point are 
shown (Figure 2a) and the proportion of antigen positive cells with 
cytoplasmic staining (for 5 images each of N= 15-70 cells/63x 
image) are shown for each time point in Figure 2b. At all time 
points, cells infected with HK-wt virus demonstrated a predom- 
inandy nuclear NSl protein localization (Figure 2a), with an 
overall average of 12.4% of positively stained cells that demon- 
strated NSl localized to the cytoplasm (Figure 2b). Ho\vc\^(t, both 
PR8 and WSN viruses had a different NSl localization phenot^'pe, 
where each induced a significandy greater percentage of infected 
cells with NSl detected in die cytoplasm at 4 h than HK-wt (33% 
and 28% versus 9.3%, respectively (p<0.01)) and again at 8 h 
(31.7% (p<0.05) and 84.3% (p<0.0()01) versus 8.4%,, respectively) 
but that were not significandy different from HK-wt at later time 
points (Figure 2b). 

NSl mutants M106I+L98S and D125G, like die PR8 and WSN 
viruses, also demonstrated an increase in cells positive for NSl in 
the cytoplasm, however the kinetics of NSl localization were 
delayed. Thus unlike PR8 and WSN, die rHK NSl mutants 
showed a temporally j)r()gTessi\'c increase in cytoplasmic positive 
NSl cells, with the greatest percxntage at 24 hpi, where 85%) or 
65% of NSl antigen positive cells stained for cytoplasmic NSl 
respectively (p<0.00001 and p<0.0001, relative to HK-wt) 
(Figure 2b). These results indicated that the mouse-adaptive 
M106I+L98S and D125G NSl mutations increased the preva- 
lence of infected cells positive for cytoplasmic NS 1 protein by a 
mechanism that increased with time following infection. 



Adaptive NSl Mutants Increased the Relative 
Cytoplasmic to Nuclear Content of the NSl Protein in 
Infected Mouse Cells 

We next quantified the effect of the adaptive NS 1 mutations on 
altered subcellular distribution of the NSl protein in infected 
mouse cells by Western blot of cellular fractions. Ml cells were 
infected with HK-wt and NS 1 mutant viruses, as well as prototype 
PR8 and WSN viruses using a MOI = 2. Following 16 hpi, die 
cells were lysed with 1 % NP-40 to acquire the whole cell fraction, 
or differentially centrifuged to obtain nuclear and cytoplasmic 
cellular fractions (n = 2-3), as per published protocols [23]. 
Figure 3a shows representative Western blots of cells infected 
with each virus (whole cell lysate (wc), nuclear fraction (n), and 
cytosolic (c) fraction) to detect NSl protein as well as cytoplasmic 
(tubulin) and nuclear (histone H3) reference standards for 1 of 2 or 
3 representative experiments. Densitometry was performed to 
quantify relative NSl abundance in each cellular fraction, which 
was normalized to the abundance of the tubulin or histone H3 
internal standards, respectively. We then determined the ratio of 
cytoplasmic to nuclear NS 1 protein abundance (c/ n NS 1 index) by 
normalizing input levels to tubulin and histone H3 band intensities 
of whole cell fractions, respectively (Figure 3b; detailed in 
methods). 

The HK-wt NSl protein was found to be predominantiy 
nuclear in mouse Ml cells (0.27 c/n NSl index) (Figure 3b), which 
was in agreement with the confocal immunostaining of mouse cells 
infected with the HK-wt virus where cytoplasmic NSl was 
detectable in the minority of antigen positive cells relative to strong 
nuclear staining (Figures 1—2). However, all NSl mutants, as well 
as the prototype mouse-adapted strains WSN and PR8, resulted in 
a trend of >2-fold increase in the c/n NSl index compared to the 
HK-wt virus (p= 1.12x10^, two-tailed t-test; Figures 3b). While 
the greatest proportion of NSl localized to the cytoplasm was 
demonstrated by both mutations at aa site 106 (Ml 06V (not 
statistically significant) and Ml 061 (p<0.05)), cytoplasmic NSl was 
reduced when expressed as double mutations, in combination with 
Ml 241 and L98S, respectively (Figure 3b). These data indicated 
that an increased proportion of NSl localized to the cytoplasm is a 
general phenotype among NS 1 mouse-adapted mutants of human 
lAV. It is possible that the low proportion of human HK-wt virus 
NS 1 protein in the cytoplasm of mouse cells may be a functional 
defect that is under selective pressure. 

To control for equivalent levels of virus infection, we also 
probed for the viral nucleoprotein (NP) and matrix protein (Ml), 
which corroborated previously published results [21], showing 
greater protein synthesis for all NSl mutants(Figure S3a). 

Altered Human Influenza virus NSl Localization Induced 
by Mouse-adapted Mutations is a Host Dependent 
Phenotype 

As an altered cellular distribution of the NSl protein in infected 
mouse cells was observed for all MA NSl mutants, we next asked 
whether this phenotype was host dependent. To this end, we 
assessed NSl locahzation in infected human A549 cells by 
immunohistochemistry and confocal microscopy for NSl mutants 
D125G and M106I-I-L98S, as well as the HK-wt, PR8 and WSN 
viruses. The cells infected with HK-wt displayed NSl staining in 
both the nucleus and cytoplasm of A549 cells. This phenotype was 
also observed for cells infected with the aforementioned NSl 
mutants as well as WSN (Figure 4a). However, cells infected with 
PR8 stained for NSl primarily in the nucleus in human A549 cells 
(Figure 4a). 
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DAPI NS1 MERGE DAPI NS1 MERGE 




Figure 1. MA NS1 mutants alter subcellular localization of the NS1 protein in infected mouse cells. Mouse IVI1 cells were infected at 
IVIOI = 3 with rHK NS mutant viruses or the HK-wt virus, or mock infected with PBS. Following 1 6 hpi, cells were fixed and stained using a polyclonal 
anti-NSI primary antibody and Cy3-conjugated secondary antibody as well as DAPI to localize the nucleus. Representative images are shown (taken 
using a 63 X oil immersion objective). 
doi:1 0.1 371 /journal.pone.0084673.g001 



We then proceeded to quantitatively determine whether the 
distribution of the NS 1 protein in infected A549 cells was altered 
by mouse-adapted NSl mutations by cell fractionation and 
Western blot (Figures 4b, S3b). The HK-wt NSl protein possessed 
a higher c/n NSl index in infected human cells compared to 
infected mouse cells (0.64 versus 0.27; Figures 3b, 4c), and the 
majority of NSl mutants (all but L98S and M106I+L98S) 
exhibited a similar or decreased proportion of cytoplasmic NS 1 



compared to HK-wt (Figure 4c). These observations suggest that 
altered NSl localization is predominantly a host-dependent 
phenotype associated with lAV host adaptation, with the 
exception of the NSl mutants possessing a mutation at L98S, 
which induced a higher proportion of NSl localized to the 
cytoplasm in both cell lines (Figures 3 and 4). 
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Figure 2. MA NS1 mutants M106I+L98S and D125G temporally increase cytoplasmic NS1 localization in infected mouse cells. Mouse 
Ml cells were infected at M0I = 3 with the rHK NS mutant viruses, the HK-wt virus, or prototype mouse-adapted viruses A/PR/8/34 (H1N1) (PR8) or A/ 
WSN/33 (H1N1) (WSN). Following 4, 8, 16, and 24 hpi, cells were fixed and stained using a polyclonal anti-NSI antibody and Cy3-conjugated 
secondary as well as DAPI to localize the nucleus. (A) Representative images are shown, taken at 63 x using oil immersion. (B) Data represent the 
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average percentage of cells positive for the NS1 antigen detected in the cytoplasm ± SD (analysis of n = 5 randomly selected images) for 4-24 hpi 
(*p<0.05, ***p<0.001, ****p<0.0001, *****p<0.00001; two-tailed student's t-test compared to HK-wt values). 
doi:1 0.1 371 /journal.pone.0084673.g002 



Increased Distribution of NSl to the Cytoplasm is 
Associated with a Host-dependent Enhancement of viral 
Fitness 

We have previously shown that MA NS 1 mutations induced an 
average of a 21 -fold increase in viral replication in mouse Ml cells 
compared to the HK-wt virus [21], and now report that these 
mutants tend to increase distribution of NSl to the cytoplasm in 
mouse cells (Figure 3). We next explored whether the host- 
dependent alteration of NS 1 localization associated with MA NS 1 
mutations correlated with viral replicative fitness by assessing virus 
replication in human A549 cells (n = 3). As the HK-wt virus was a 
human isolate, and presumably adapted to the human host, we 
hypothesized the virus would replicate to a higher level in human 
cells than in mouse cells. In respect to previously published virus 
growth curves in infected mouse cells, the HK-wt virus maximum 
yield was >1500 fold greater in human cells than in mouse cells 
(5.9x10*^ pfu/mL versus 3.8x10^ pfu/mL [21], (Figure 5). This 
difference in host range for the HK-wt virus was consistent with 
diminished NP protein expression in human cells (Figure S3c). 

The majority of viruses expressing NSl mutations sc'kxted upon 
mouse adaptation were unable to increase virus replication in 
human cells, with 9 of the 1 2 mutants replicating to levels S HK- 
wt (Figure 5). While NSl mutants Ml 061, and to a lesser extent, 
F103L and M106V, did in fact increase virus replication 
throughout the time course (p<0.001, p<0.05, and p<0.05, 
respectively), the average level of enhancement among NSl 
mutants amounted to only 1.5 fold; Figure 5). These results suggest 
that replicative fitness acquired upon adaptation to a novel host is 
largely a host-dependent phenotype for the NSl mutations tested. 

Altered Host Gene Expression in Mouse Cells following 
lAV Infection due to Mouse Adapted NSl Mutations 

The NS 1 protein has been shown to bind two cellular proteins 
involved in the post-transcriptional polyadenylation of host gene 
transcripts; nuclear poly-A binding protein II (PABPII), and 
cleavage and polyadenylation specificity factor 30, (GPSF30). The 
ability of NSl to bind CPSF30 has been shown by transfection to 
positively correlate with reduced gene expression of target genes 
including IFN P [2,26]. We have previously shown that MA NSl 
mutations induced significant changes in binding affmity for 
CPSF30, with 8 mutants losing or reducing binding capacity [21]. 
However, there was no correlation between the ability to bind 
CPSF30 and the ability to shut off mouse IFN-P production in vivo 
[21]. We next wanted to determine whether the ability of the NSl 
protein to bind CPSF30 was associated with a more general 
regulation of gene expression in mouse cells. To this end we 
performed microarray analyses using the mouse Affymetrix 
Genechip Mouse Exon 1.0 ST Array on RNA extracted from 
mouse cells 8 hours following infection with HK-wt and NSl 
mutants at a MOT of 2. Three independent infections were 
analyzed in parallel with comparison of genes that were 
significandy up or down regulated by ^ 1 or S 1 log 2 fold (up 
or down regulated by >2.0 fold; p<0.05 by ANOVA (n = 3)) 
relative to mock-infected Ml cells. On determining the host genes 
that were differentially regulated relative to mock infected cells, we 
observed that the viruses possessed one of 2 phenotypes, which 
had either low or high effects on host gene expression (Figure 6a; 
p<0.05, number of genes 2 fold up- or down- regulated in 
expression compared to mock-infected cells). The low gene- 



regulating (LGR) group consisted of HK-wt and 5 of the 12 NSl 
mutants (D125G, M106I, V18()A, V226I, and R227K), where the 
number of genes up regulated exceed the number down regulated 
(pS0.05), (average of 53 genes up regulated and 8.5 genes down 
regulated) (Figure 6a). Mouse cells infected with the 7 remaining 
MA NSl mutants (D2N, V23A, F103L, M106I+L98S, L98S, 
M106V, and M106V4-M124I) were found to activate more gene 
regulation events; compared to the LGR group of viruses, this high 
gene-regulating group (HGR) induced up-regulation and down- 
regulation of 8-fold and 93-fold more genes, respectively (p£0.05) 
(Figure 6a). 

We then performed hierarchical analyses of gene expression for 
the 5274 genes that were differentially regulated by at least one 
mutant. Gene regulation differences were illustrated upon 
histogram analysis of the HK-wt NS 1 and each mutant to show 
the distribution of up and down regulated genes for each mutant 
(Figure S4). The six members of LGR group clustered together 
according to differential gene regulation. The M106I and R227K 
NSl mutants most closely resembled HK-wt- induced host gene 
regulation, and differed from the closely linked V180A and V226I 
that were secondarily linked to the D125G mutant (Figure 6b). 
The hierarchical clustering of the HGR group showed 2 closely 
related clusters. NSl mutant V23A clustered with M106I+L98S, 
and this hierarchical grouping was linked secondarily to a second 
cluster of mutants; F103L and Ml 06V, which both previously 
were observed to induce similar host gene expression profiles in 
infected mouse cells (Figure 6b) [24]. The differentially regulated 
genes also fell into two main clusters of up and down regulated 
genes that were further divided into subgroups that were shared to 
different extents amongst the mutants (Figure 6b), indicating two 
distinct mechanisms of gene regulation for the NSl mutants. 
Comparison of the gene regulation patterns (given as a relative 
statistic of the sum of the absolute values of the log2 fold 
differentially regulated genes for each virus) and virus properties of 
virus yields in Ml cells and CPSF binding [21], show a general 
lack of relationship among these parameters (Figure 7). Although 
the HGR group tended to generate higher yields of virus in mouse 
cells than the LGR group, both groups possessed high and low 
yielding viruses. These results suggest that the effects on host gene 
expression are not simply a measure of the extent of replication or 
ability to bind CPSF but are complex parameters (Figure 7). 

Gene ontolog}' (GO) analysis was performed at Level 2 for the 
average gene expression levels of the differentially regulated genes 
in the LGR and HGR groups. Although both groups included 
genes that were up and down regulated for host genes involved in 
biological processes, cellular components, and molecular func- 
tions, they differed fundamentally in the overall effect of these 
changes. The LGR group effects were dominated by positive 
regulation of host genes whereas the HGR group was dominated 
by negative regulation of genes suggesting that they provide 
adaptive gain of function through different mechanisms of host 
interaction (Figures So, S6). 

The HGR and LGR mutations mapped to two regions of the 
NSl protein, divided at aa position 106 (Figure 7 and 8). While 
these regions do not correspond with the accepted functional 
domain partitioning for the NSl protein (RNA binding domain 
(aa 1-73) and effecter domain (aa 74—230) [1], our results do 
suggest different roles of the amino terminal versus carboxy- 
terminal domains of the NSl gene in regulating host gene 
expression. 
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Figure 3. MA NS1 mutations alter subcellular localization of the NS1 protein in infected M1 ceils. Mouse iVI1 cells were infected (M0I = 2) 
with rHK NSl mutant viruses, the HK-wt virus, or the PR8, or WSN viruses, or mock infected with PBS. At 16 hpi the cells were lysed and differentially 
centrifuged to obtain whole cell (w), nuclear (n) and cytosolic (c) fractions. Cell fractions were resuspended in SDS buffer then were separated by SDS 
PAGE electrophoresis, followed by Western blot analysis to detect the NSl protein as well as loading markers tubulin and histone H3. (A) 
Representative blot of 2-3 independent experiments, which were (B) analysed by densitometry where the amount of NSl or NS3 protein in the 
nuclear and cytoplasmic fractions of a cell infected with a given virus was normalized to histone H3 and tubulin levels of the whole cell lysate, 
respectively. Cytoplasmic distribution was calculated by dividing cytoplasmic NSl by nuclear NSl. Data represent the means ± SE (*p<0.05, 
**p<0.01; two-tailed student's t-test). 
doi:1 0.1 371/journal.pone.0084673.g003 
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Figure 4. Increased cytoplasmic distribution of NS1 is a Kiost-dependent pKienotype of MA NS1 mutant viruses. (A) Human A549 cells 
were infected at M0I = 3 with the rHK NS mutant viruses, the HK-wt virus, or the WSN or PR8 viruses. Following 16 hpi, the cells were fixed and 
stained using a polyclonal anti-NSI antibody and Cy3-conjugated secondary antibody as well as DAPI to localize the nucleus. Representative images 
are shown, taken at 63 x using oil immersion. (B) Human A549 cells were infected with rHK NS mutant viruses, the HK-wt virus, or the WSN or PR8 
viruses, or mock infected with PBS. At 16 hpi cells were lysed and differentially centrifuged to obtain whole cell (w), nuclear (n) and cytosolic (c) 
fractions. Cell fractions were resuspended in SDS buffer then were separated by SDS PAGE electrophoresis, followed by Western blot analysis to 
detect the NS1 protein as well as loading markers tubulin and histone H3. (B) Representative blot of 2-3 independent experiments, which were (C) 
analysed by densitometry where the amount of NS1 or NS3 protein in the nuclear and cytoplasmic fractions of a cell infected with a given virus was 
normalized to histone H3 and tubulin levels of the whole cell lysate, respectively. Cytoplasmic distribution was calculated by dividing cytoplasmic 
NS1 by nuclear NS1. Data represent the means ± SE (*p<0.05, two-tailed student's t-test). 
doi:1 0.1 371 /journal.pone.0084673.g004 



NS1 Mutants Enhance Genomic viral RNA Replication in 
Infected Mouse Cells 

We have previously reported the majority of NSl mutants 
significantly increase viral mRNA abundance in infected mouse 
cells [21]. We extended this analysis to determine the effects of the 
NSl mutants on genomic viral RNA (vRNA) transcription. RNA 
extracts employed for microarray analysis were subjected to real- 
time quantitative RT-PCR analysis of genomic viral RNA (vRNA) 
encoding NSl, NP and M 1 genes (n = 3 biological replicates with 2 
technical replicates of each qPCR analysis). All NSl mutants with 
the exception of V180A and V226I significantly increased the 
abundance of NSl, NP and Ml vRNA genes relative to HK-wt 



levels (pS0.05, paired t-test; Figure 9). Comparing the fold 
differences relative to HK-wt levels, relative abundance of vRNA 
for the three tested transcripts were similar for V23A, L98S, 
F103L, M106I, M106V, V180A, and V226I, but differed for the 
remaining mutants (D2N, M106I-HL98S, and R227K), that tended 
to possess higher relative abundances of the M 1 vRNA relative to 
the NP and NSl genome segments (Figure 9). These data 
demonstrated that the majority of NSl mutations selected upon 
LAV adaptation to the mouse increase both viral transcription [21] 
and genome replication (Figure 9) but to different extents and 
different patterns among the three genes tested. 
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Figure 5. MA NS1 mutations enKiance virus growth in human 
A549 ceils in vitro. Human A549 cells were infected in triplicate at an 
IVIOI of 0.02 with rHK NS mutants or HK-wt virus, and supernatant 
collected 12, 24, 48, and 72 hpi was quantified for viral litre by plaque 
assay on MDCK cells. Data represent the mean viral titre of n = 3 
replicates (each titrated in duplicate) ± SD (♦p<0.05, ♦*p<0.01, 
***p<0.001; two-tailed student's paired t-test for titers obtained from 
12 to 72 hpi). 

doi:10.1371/journal.pone.0084673.g005 
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Figure 6. Mouse-adaptive NS1 mutations regulate mouse gene 
expression by two different mechanisms. Mouse Ml cells were 
infected in triplicate at an MOI of 2 with rHK NS mutants or HK-wt virus, 
and total RNA isolated at 8 hpi was analyzed by microarray using 
GeneChip Mouse Exon 1.0 ST Array (Affymetrix, Santa Clara, CA, USA). 
Microarray gene expression data were normalized and analyzed by 
Flexarray 1.6.1. Genes were considered as up or down regulated relative 
to mock infected cells if they were <1 or >1 log2 fold different (£ or 
£2 fold differences) expression level (p<0.05; ANOVA). (A) The number 
of differentially regulated host genes that were significantly up or down 
regulated £2 or 2:2 fold at the p<0.05 by ANOVA is plotted for each 
mutant. The mutants formed two groups with either "low gene 
regulation" (LGR) or "high t gene regulation" (HGR) phenotypes. (B) 
Heat map of differentially regulated genes in mouse cells infected with 
HK-wt or rHK NS1 mutant viruses relative to mock infected Ml cells. 
Genes (total = 5274) were included for hierarchical clustering analysis 
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among mutants if they were differentially regulated (£2~^or >:2^ fold 
differences) and significantly different from mock infected cells for one 
or more of the mutants, and gene regulation signatures were also 
analyzed for hierarchical clustering among viruses. The scale depicts up 
(red) and down (blue) regulated genes according to the log2 scale 
shown; equal expression is indicated in white (log 2°= 1). 
doi:l 0.1 371/journal.pone.0084673.g006 



Phenotypic Trends of Mouse Adapted NSl Proteins 

Thus far we have shown several novel adaptive phenotypes 
associated with NSl mutations selected on mouse adaptation, 
including altered NSl protein localization, enhanced vRNA 
production, and regulation of host gene expression. Previous 
reports have indicated the ability to bind CPSF30 correlated with 
nuclear localization of the NSl protein [12] as well as ability to 
shut down host gene expression by inhibiting pre-mRNA 3' 
processing [27]. We tabulated the known phenotypic properties of 




Figure 7. NS1 mutant CPSF30 binding affinity and ability to enhance viral replication in mouse cells does not correlate with 
regulation of mouse gene expression. Top panel: maximum virus yield in mouse Ml cells infected using a MOI of 0.02 relative to HK-wt virus 
yield (relative growth, top panel; values from [21]); middle panel: binding affinity to the CPSF30 F2F3-FLAG protein relative to HK-wt NSl (normalized 
CPSF30 binding, middle panel; values from [21]), bottom panel: the cumulative effect on host gene expression of both up and down regulated genes 
as the sum of the absolute values of the log2 expression levels of host genes for each NSl mutant. 
doi:1 0.1 371/journal.pone.0084673.g007 
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Figure 8. Protein map of NS1 indicating the location of mutations selected upon HK virus mouse adaptation. Linear map of NS1 
protein indicating the location of NS1 mutations selected upon adaption of A/Hong Kong/1/1968 (H3N2) to the mouse, as depicted by black lines. 
NS1 nuclear localization sites (NLS1 and 2), nuclear export signal (NES), and its adjacent negative regulator of NES are indicated on the protein map as 
well as the amino and carboxyl protein regions associated with sites of mutations inducing high or low host gene regulation phenotypes (HGR or 
LGR, respectively). 

doi:1 0.1 371 /journal.pone.0084673.g008 



the NS 1 mutant panel to show that ability to bind the host factor 
CPSF did not correlate with intracellular localization of the NS 1 
protein (Figure 10). We therefore conclude additional mechanisms 
beyond CPSF binding are involved in host gene regulation by the 
NSl gene. 

Discussion 

We have previously shown that the NSl mutations selected on 
lAV evolution in the mouse are adaptive, possessing multiple gain- 
of-function phenotypes including enhanced viral replication in vivo 
in the mouse lung and in vitro in mouse cells, enhanced viral gene 
expression at both the level of viral mRNA and viral protein in 
mouse cells, and ability to suppress IFN-P production in the mouse 
lung [21]. Moreover, observed reductions in IFN-P production 
in vivo did not correlate with increased binding to the cellular post- 
transcriptional processing factor CPSF30 [21] (Figure 10). In this 
paper, we expand the scope of phenotypes associated with MA 
NS 1 mutations and demonstrate that host adaptation of lAV is 
associated with increased distribution of the NS 1 protein into the 
cytoplasm of infected cells (Figures 3 and 4), and that this 
phenotype is host-dependent and thus restricted to the adapted 
host (mouse), and not the original host (human) (Figure 10). We 
also demonstrated that NSl mutations selected upon mouse 
adaptation altered host gene expression profiles by one of two 
distinct but unknown mechanisms. Notably, gene ontology 
analysis showed that while the LGR phenotype was associated 
with positive regulation of genes involved in host processes, these 



host processes were negatively regulated for the HGR phenotype 
(Figures 6, S5, and S6). 

Altered NSl Cellular Distribution is Associated with 
Mouse Adaptive Phenotypes of NSl Mutations 

The NS 1 protein of the human prototype clinical isolate, HK- 
wt, was abundant in both the cytoplasm and nucleus of infected 
human A549 cells (Figure 4c) suggesting that this is an optimal 
distribution for a virus adapted to its host. Consistent with a 
defective ability to partition into the cytoplasm in the novel mouse 
host cell, HK-wt was predominately associated with the nucleus 
(Figure 3b). However all 12 mutant NSl genes expressed in 
recombinant viruses in mouse cells showed a tendency to an 
increased proportion of NS 1 protein in the cytoplasm versus the 
nucleus with statistically significant increases for the F103L, 
Ml 061, M106I+L98S, M106V+M124I, and D125G mutants 
(Figwe 3). Prototype mouse-adapted laboratory strains PR8 and 
WSN expressed NSl proteins that accumulated in the cytoplasm 
at levels & nuclear NS 1 levels for both mouse and human infected 
cells. This is possibly due to evolutionary changes selected on 
laboratory culturing as they both possess adaptive NS 1 mutations 
23A and 2261 [21] and in addition PR8 possesses mutations 103S 
and 1061 (both of which increase virulence) [24]. It is thus likely 
that host-specific cellular distribution of the NSl protein 
contributed to the adaptive phenotypes associated with NSl 
mutations. MA NSl mutations enhance viral gene expression, 
which has previously been quantified at both the protein and 
mRNA, and now the vRNA, level in infected mouse cells (Figures 9 
and 10) [21]. Increased NSl gene expression accounts for greater 
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Figure 9. MA NS1 mutations increase genomic viral RNA syntliesis in infected mouse cells. Mouse IVI1 cells were infected in triplicate at 
an MOI of 2 with rHK NS mutant viruses or the HK-wt virus, and total RNA isolated at 8 hpi was reverse transcribed using primers specific for genomic 
viral RNA segments. Real-time PCR (qPCR) was performed to quantify NP, Ml and NS1 genomic RNA levels. Results were normalized by p-actin levels, 
and presented as values relative to HK-wt genomic RNA levels. Data represent the means ± SD (two tailed student's paired t-test) for NP, NS1 and Ml 
transcript relative levels (indicated by bracket) or two-tailed student's t-test for individual mRNA samples (n = 3) relative values (i p<0.05; ii p<0.01; iii 
p<0.001; iv p<0.0001; v p<0.00001). 
doi:1 0.1 371 /journal.pone.0084673.g009 



cellular abundance of NSl proteins; however the altered 
distribution to favour cytoplasmic localization suggests potential 
differences in nucleo-cytoplasmic transport or binding to host 
proteins, the extent to which remains unknown. 

Following its synthesis by host translation machinery, nascent 
NSl is found in the host cell cytoplasm, then rapidly localizes to 
the nucleus [5]. Most NSl proteins, including the HK NSl, 
possess two NLS, NLSl at aa 34-38 and NLS2 at aa 203-237 
[28], which by interaction with importin-ot [29], function in 
driving nuclear localization of NS 1 . The NS 1 protein also contains 
one characterized NES (aa 138-147 FDRLETLILL) that may 
contribute to the cytoplasmic localization of the NSl protein 
reported at later times following infection [6]. Regulation of 
nuclear export appears to be complex, as the adjacent region to 
the NES (aa 148-161) inhibits NES activity in transfected cells but 
not in infected cells [7]. Cellular localization of NSl protein thus is 
a culmination of multiple counteracting activities such as NLS 
(1-1-2) versus NES (+regulator activities) that shuttle NSl between 
the two ceUular compartments, and this phenomena is also 
influenced by NS 1 protein mocUfications or interactions with other 
viral proteins as described for several other viral proteins [30,3 1] . 

In further support of our findings, NS 1 cytoplasmic localization 
has also been mapped to aa 103 and 106 residues of the PR8 virus 
in MDCK cells [12], and the MIOII residue mutation in 2003 
HPAI H5N1 in infected A549 cells, that corresponds to the 106 
position in the HK-wt virus [32]. NSl aa residues 103L and 1061 
have also been linked to cytoplasmic localization of the H5N1 
highly pathogenic avian influenza (HPAI) isolate A/HK/156/ 
1997 NSl protein in mouse cells, when expressed on the PR8 viral 
backbone [24]. With respect to these findings, it is possible the 
NS 1 protein contains an additional domain that influences nuclear 
export or cytoplasmic retention spanning from aa 98 to 106. 



Additional experiments on the localization of NSl and NS3 
proteins by site directed mutagenesis of aa 98-106 or transference 
of this motif (LSEDWFMLM) to other proteins, is required to 
further define the functional roles of component aa in nuclear 
export and cytoplasmic localization. Nonetheless, the observation 
of adaptive mutations in the NSl gene that affect subcellular 
distribution is consistent with the common theme of adaptive 
evolution of nuclear localization signals in RNA polymerase 
subunits that overcome defects in interacting with cellular 
trafficking proteins [15,18,33,34]. 

Functional Roles of Enhanced Expression and Subcellular 
Localization of the NSl Protein 

We have shown that NSl adaptive mutations increase viral 
transcription and viral gene expression at the level of protein 
synthesis (Figures 3 and 10; [21]). Because NSl functions in 
enhancing translation of viral mRNA [35,36], mutants with 
increased accumulation of NS 1 in the cytoplasm of infected mouse 
cells may further promote protein synthesis. Given the NSl 
protein is more abundant in the cytoplasm of infected mouse cells 
for the majority of NSl mutants compared to HK-wt NSl, 
cytoplasmic NS 1 can more effectively promote initiation of viral 
protein synthesis by binding to viral mRNA and recruiting it to 
cellular translational factors or alternatively that NSl mutants 
favourably affect binding affmity for translational factors such as 
eIF4GI and PABPI [35,36] . Binding to cytoplasmic factors such as 
these could also contribute to sequestration of the NS 1 protein into 
the cytoplasm. Nonetheless, altered host-dependent cellular 
localization of NSl mutants could contribute positively to the 
enhancement of viral gene expression. 
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Figure 10. Adaptive properties of rHK viruses expressing NS1 mutations selected upon mouse adaptation. Heat map illustrating 
adaptive phenotypes (Figures 3-5, 7, 9, and [21]); relative to HK-wt values, where 1 is indicated as white* and the red scale depicts phenotypes 
greater than HK-wt (increasing shades of red from 1 to 10 (with saturated red for a 10), and shades of blue scale from 1 to 0 depicts phenotypes lesser 
than HK-wt values, where data was calculated as a fold change relative to HK-wt value. For virus yield values, maximum virus yield obtained in the 
growth curve was used to calculate fold change relative to HK-wt virus maximum virus yield. *With the exception of % mortality in the mouse where 
0% mortality is indicated as blue and any mutation inducing mortality in the mouse indicated as red. 
doi:1 0.1 371 /journal.pone.0084673.g01 0 



Interferon antagonism was another NS 1 property enhanced by 
the majority of NSl mutants (Figure 10, [21]), and likewise this 
adaptive phenotype may be partially attributed to altered cellular 
localization of the NSl protein. In the cytoplasm of infected 
human cells, NSl inhibits IFN production at the transcriptional 
level of RIG-I signalling, by directly binding RIG-I and multiple 
RIG-I activators including E3 hgase, TRIM25 [37], and dsRNA 
[38-40]. RNA binding as well as RIG-I and the E3 hgase Riplet 



binding [41] inhibit signaling and IFN 
a recent report by Dankar et al. (201 
1061 of the A/HK/156/1997 (H5N1) 
increase NS 1 localization to the cytopl 
and also to mediate increased binding 
with decreased IFN-fi production in 
Thus the increased distribution of NS ] 



activation in mouse cells. In 
3), NSl residues 103L and 
HPAI strain were found to 
asm of infected mouse cells, 
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infected mouse cells [24]. 
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for MA NS 1 mutants may also contribute to inhibition of type I 
IFN production. 

In addition, increased cytoplasmic NSl may also contribute to 
an enhanced ability of NSl to circumvent the IFN response. 
Cytoplasmic NSl binds dsRNA, which sequesters it from 
activation of PKR and 2'5'-OAS [10,42^4]. NSl also inhibits 
dsRNA-independent PACT-mediated PKR activation [42]. 
Moreover, cytoplasmic NSl is responsible for blocking STAT-1 
signalling induced by type I IFN, which results in suppression of 
IFN-stimulated gene transcription [2] . We have previously shown 
that the majority of MA NSl mutants shut down IFN production 
in vivo in infected mouse lungs (Figure 1 0) and suppress the effects 
of IFN on viral replication and viral protein synthesis in mouse 
cells [21]. However additional experiments including quantifica- 
tion of NSl cellular distribution at earlier times following infection 
are required to support the hypothesis that increased distribution 
of NSl to the cytoplasm increases NSl IFN antagonistic 
properties. 

Alternative Mechanisms for NSl- mediated Host Gene 
Induction and Suppression 

Gene ontology analysis of mouse genes regulated upon 
influenza virus expression indicates two general mechanisms of 
regulation of host gene expression that involve alternatively down- 
or up regulation of host gene expression that mapped to the amino 
or the carboxyl- regions spanning aa 2 to 106 and 106 to 227, 
respectively (Figures 6, 7, 8, 10 and S5-S6). This suggests that NSl 
may be interacting with transcriptional regulators or factors in the 
nucleus and or cytoplasm, and that mutations selected upon 
mouse adaptation may be altering such interactions. Human 
H3N2 NSl genes isolated since 1969 have been shown to possess a 
histone-H3 mimic sequence ARSK at position 226-229 that 
functions to competitively inhibit activation of host anti-viral 
mediator genes to enhance virus replication [45]. The 226-229 
VRSK sequence in HK-wt NSl was mutated in the V226I and 
R227K mutants raising the possibility that these mutations affect 
gene interactions through modification of this system; confirma- 
tion of which will require further analysis. 

The NS 1 protein has been reported to block host gene expression 
at the level of mRNA production by interaction with numerous host 
factors. In addition to binding CPSF30, NSl also binds the polyA 
binding protein nuclear 1 (PABPN 1 ) [46] , as well as the host mRNA 
nucleocytoplasmic export machinery, NXFl/TAP [47], which 
leads to the inhibition and thus down regulation of nuclear export of 
processed mRNA. Further assays of protc-in and RNA interactions 
between NSl mutants and other nuclear cellular binding factors 
would be necessary to determine whether altered affinity to host 
gene regulation factors is involved in the mechanism by which the 
NSl protein alters host gene expression. Given the set of NSl 
mutants tended toward one of two distinct gene regulating 
phenotypes, this suggests more than one NSl gene-regulating 
mechanism is selected upon influenza virus adaptation to a novel 
host which may be further dependent on other adaptive mutations 
that are selected in other viral genes [1-4,14]. 

In summary we have demonstrated a deficiency in cytoplasmic 
NS 1 as a host barrier to infection with a human influenza virus in 
mouse cells that was o\XT('()m<; by adapti\'e mutations. Although 
there have been no previous reports on the roles of host adaptation 
in modulating NSl subcellular localization, these findings dem- 
onstrate that host adaptation alters the ability of the NS 1 protein 
to localize to the cytoplasm in infected ceUs, and this altered 
ceUular distribution was associated with other adaptive phenotypes 
such as enhanced IFN antagonism and viral gene expression. 
Moreover, we have shown NSl mutations selected upon mouse 



adaptation alter host gene regulation by different, but unknown, 
mechanisms irrespective of the ability to bind CPSF30. 

Supporting Information 

Figure SI Infection with NSl mutants results in detect- 
able levels of NP production. Mouse M 1 cells were infected at 
MOI = 3 with rHK NS mutant or rHK-wt viruses as indicated. 
FoUowing 16 hpi, cells were fixed and stained using a polyclonal 
anti-NSl antibody and Cy3-conjugated secondary as well as DAPI 
to locahze the nucleus. Representative images are shown, taken at 
63 X using oil immersion. 
(TIF) 

Figure S2 Quantification of virus-infected mouse cells 
with detectable NSl protein in the cytoplasm. Mouse Ml 
ceUs were infected as described in figure 1. Data represent the 
average percentage of ceUs positive for the NSl antigen detected in 
the cytoplasm ± SD (analysis of n = 5 randomly selected images) at 
16 hpi (*p<0.05, ***p<0.001; two-tailed student's t-test compared 
to HK-wt values). 
(TIFF) 

Figure S3 Viral protein expression in mouse and 
human cells for rHK NSl mutant or HK-wt viruses. 

Cellular fractions of infected mouse cells (A; same samples as 
shown in figure 3) or infected human whole cell lysates (B; same 
samples as shown in figure 4) were separated by SDS page then 
probed with anti-NP and anti-Ml antibodies. (C) HK-wt virus 
produces less viral protein in infected mouse ceUs than in infected 
human cells (same samples as in figures 3 and 4). WC: whole cell 
lysate; N: nuclear fraction; C: cytoplasmic fraction; NP: nucleo- 
protein; Ml: matrix protein 1. 
(TIFF) 

Figure S4 Histogram analysis of the HK-wt NSl and 
each mutant to show the distribution of up and down 
regulated genes. The relative distribution of number of genes 
and their log2 regulation levels relative to mock infected Ml cells 

are shown for the genes analyzed by hierarchical clustering in 
Figure 6. The up and down regulated genes are show in shades of 
red and blue respectively according to the scale shown with values 
of 1 shown in white (2" = 1). 
(TIFF) 

Figure S5 Gene ontology (GO) analysis was performed 
at Level 2 for the average gene expression levels of the 
differentially regulated genes in the low host gene 
regulating (LGR) group. The average gene expression levels 
among the set of 5,274 differentially regulated genes of the 
member of the LGR group of mutations (D125G, Ml 061, V180A, 
V226I, and R227K) were characterized for their up or down 
regulating effects for host genes involved in biological processes, 
ceUular components, and molecular functions. The data show the 
number of genes that are either up or down regulated for each of 
the functional areas indicated. The LGR group effects were 
dominated by positive regulation of host genes. 
(TIFF) 

Figure S6 Gene ontology (GO) analysis was performed 
at Level 2 for the average gene expression levels of the 
differentially regulated genes in the high host gene 
regulating (HGR) group. The approach described for Fig S2 
was done for the HGR group of mutants, p2N, V23A, F103L, 
M106I+L98S. L98S, M106V, and M106V+M124I). The HGR 
group was dominated by negative regulation of genes. 
(TIFF) 
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